Discrete and mesoscopic regimes of finite-size wave turbulence 
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Bounding volume results in discreteness of eigenmodes in wave systems. This leads to a depletion 
or complete loss of wave resonances (three-wave, four-wave, etc.), which has a strong effect on 
Wave Turbulence, (WT) i.e. on the statistical behavior of broadband sets of weakly nonlinear 
waves. This paper describes three different regimes of WT realizable for different levels of the wave 
excitations: Discrete, mesoscopic and kinetic WT. Discrete WT comprises chaotic dynamics of 
interacting wave "clusters" consisting of discrete (often finite) number of connected resonant wave 
triads (or quarters). Kinetic WT refers to the infinite-box theory, described by well-known wave- 
kinetic equations. Mesoscopic WT is a regime in which either the discrete and the kinetic evolutions 
alternate, or when none of these two types is purely realized. We argue that in mesoscopic systems 
the wave spectrum experiences a sandpile behavior. Importantly, the mesoscopic regime is realized 
for a broad range of wave amplitudes which typically spans over several orders on magnitude, and 
not just for a particular intermediate level. 



INTRODUCTION 

Dispersive waves play a crucial role in a vast range of 
physical applications, from quantum to classical systems, 
from microscopic to astrophysical scales. For example, 
Kelvin waves propagating on quantized vortex lines pro- 
vide an essential mechanism of turbulent energy cascades 
in quantum turbulence in cryogenic Helium [1-6]; wa- 
ter waves aid the momentum and energy transfers from 
wind to ocean [7]; internal waves on density stratifica- 
tions and inertial waves due to rotation are important 
in turbulence behavior and mixing in planetary atmo- 
spheres and oceans [8-10]; planetary Rossby waves are 
important for the weather and climate evolutions [11]; 
and Alfven waves are ubiquitous in turbulence of solar 
wind and interstellar medium [12-17]. More often than 
not, nonlinear interaction of different wave modes is im- 
portant in these and other applications, and there has 
been a significant amount of work done in the past to 
describe evolution of such interacting wave systems. If 
the number of excited modes is large, they experience 
random evolutions, which must be described by a sta- 
tistical theory. Weak Wave Turbulence refers to such a 
statistical theory for weakly nonlinear dispersive waves in 
unbounded domains [18]. This approach was initiated by 
Peierls in 1929 to describe phonons in anharmonic crys- 
tals [19], and it was reinvigorated in 1960's in plasma 
physics [22-24] and in the theory of water waves [20, 21]. 
By now, it has been applied to description of a great vari- 
ety of physical phenomena, from synoptic Rossby waves 
[25-28] to magneto- hydrodynamic turbulence [14-16], to 
acoustic waves [29], to waves in stratified [8, 9] and ro- 
tating fluids [10], and many other physical wave systems. 

On the other hand, it has become increasingly clear 
that in the majority laboratory experiments and numer- 
ical simulations of nonlinear dispersive wave systems the 
discreteness of the wave-number space due to a finite size 
is a crucially important factor which causes the system 
behave differently from the predictions of the classical 
theory of wave turbulence based on the continuous (in- 



finite domain) limit [30-38]. Moreover, similar behavior 
often occurs in nature when waves are bounded, e.g. for 
planetary Rossby waves bounded by the finite planet ra- 
dius [39]. 

Description of transition from regular to random 
regimes and characterization of the intermediate states 
where both regular and random wave motions are present 
and mutually inter-connected, is an intriguing and chal- 
lenging problem. Such intermediate states where the 
number of waves is big and yet the discreteness of the 
wavenumber space still remains important are called Dis- 
crete and Mesoscopic Wave Turbulence. 



I. WEAKLY INTERACTING WAVES 

A. Normal modes of linearized problem 

An evolution equation is called dispersive if its linear 
part has wave-like solutions ip{r,t), that depend on the 
coordinate in the d-dimensional physical space, r € R'' 
and time t as follows: 



= \ Ak I cos(fc ■ r — ujt + Lp) 



(la) 



where "c.c." means "complex conjugate" . Here Ak — 
\Ak \ exp{iip) g C is a constant wave amplitude, cp is the 
wave phase, k G M."^ is a wavevector and wave frequency 
u} = io{k) e K is such that \d'^w / dkidkj\ ^ 0, where 
ki and kj are components of k. Physically, the latter 
condition means that wave packets with different wave- 
numbers propagate at different speeds, so that localized 
initial data would disperse (spread) in space. 

In bounded systems, the set of normal wave modes 
becomes discrete. For waves in a periodic d-dimensional 
cube with side L, the normal modes are given by (1) with 
a discrete set of wavenumbers k — 2ii where I G 
For different boundary conditions, normal modes of the 
linearized problem may differ from the propagating plane 
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waves (la). For instance, zero boundary conditions in 
a rectangular box typically (but not always!) lead to 
standing waves, 

tp{r,t) = |Afc| sin(A; • r + fsp) sm{ujt + ipt) , (lb) 

where (psp and (pt are the space and time phases corre- 
spondingly. A more complex form of the normal mode 

is given by ocean planetary motions in a rectangular do- 
main [0,^2,] X [0,Ly] with zero boundary conditions, see 
e.g. [50]: 



'fit 



^3\r,t) = \Ak\ sm(^7r— j sm(^7r— j sm(^— 

X y 

(Ic) 

where m, ii fc N are integers and ui = 
(3/[2ny/{m/Lx:)'^ + (n/Ly)^] with a constant /3 called 
Rossby number. 



B. Equation of motion 

A rather general class of non-dissipative nonlinear 
waves can be described within the framework of the clas- 
sical Hamiltonian approach. This means that after a 
proper change of variables the motion equation in nat- 
ural variables (fluid velocity, electrical fleld, density vari- 
ations, Cite.) can be presented in the universal form of 
canonical Hamiltonian equations for canonical variables 
b{r,t), b*{r,t), which characterize the wave amplitudes. 
Here " * " denotes complex conjugation. The Hamilto- 
nian equations for the space-homogeneous systems are 
most conveniently written in Fourier space, because it 
is a natural space for describing the wave solutions. In- 
troducing the Fourier transform of b{r,t) and calling it 
Qk = a{k,t), the Hamiltonian equation can be written as 
follows [18], 



.dttk 
I—— 
dt 



dn 

dal 



(2) 



Hamiltonian "H = 'H{ak,a'^} is usually (but not neces- 
sarily) is the energy of the wave system, expressed in the 
terms of the canonical variables Ofc, for all allowed by 
the boundary conditions wave vectors k. In the simplest 
case of a periodical box fe = 2-1^1/ L, with wavenumber 
i e Z"^ and L being the box size and d is space dimen- 
sion. 

For the waves of small amplitudes (for example, when 

the elevation of the gravity waves on the water surface 
is smaller then the wavelength) the Hamiltonian can be 
expanded in powers and a?: 



H — 'H2 + ^int ) 



(3a) 
(3b) 



where 'Hj is a term proportional to product of j ampli- 
tudes afc and the interaction Hamiltonian "Hint describes 
the wave coupling, as explained below. We omitted here 



the independent of and a^, part of the Hamiltonian 
T-Lq, because it does not contribute to the motion equa- 
tion (2). In this paper we consider only waves exited 
about steady equilibrium states, i.e. if absent initially, 
the waves must remain absent for all time, afc = a)^ = 0. 
Thus, the linear Hamiltonian is zero, "Hi = 0. 

Expansion (3b) utilizes the smallness of the wave am- 
plitudes, therefore, generally speaking. 



(4a) 



In particular cases, due to specific symmetries of a prob- 
lem, the odd expansion terms vanish (i.g. for spin waves 
in magnetics with exchange interactions, Kelvin waves 
on quantum vortex lines). In these cases, instead of (4a) 
one has: 



7^3 — T~L^ — Hi — • • • — , 

^4 > 7^6 > ^8 > . . . . 



(4b) 
(4c) 



Three- wave interactions often dominate in wave systems 
with small nonlinearity, e.g. for Rossby waves in the At- 
mosphere and Ocean, capillary waves on the water sur- 
face, drift waves in plasmas, etc. On the other hand, if 
"Ha = 0, or if three-wave resonances are forbidden (in 
the sense that will be clarified below) the leading non- 
linear processes may be four- wave interactions. Further, 
there are examples of systems where the four-wave in- 
teraction is absent and the leading nonlinear process is 
five-wave, e.g. for one-dimensional gravity water waves 
[42-44], or even six-order, e.g. for Kelvin waves on quan- 
tum vortex lines [3, 45] . However, such higher-order wave 
systems are rather rare and, therefore, in this paper we 
will discuss three- and four- wave interactions only, which 
describe most of weakly interacting waves. 



C. Non-interacting waves 

The first physically meaningful expansion term, 
quadratic Hamiltonian 



■H2 = y^tjfclafcl 



(5a) 



n=l 



according to Eq. (2) produces a linear equation of motion. 



dak 
~dt 



UJkdk ■ 



(5b) 



and thus describes noninteracting waves with the disper- 
sion relation Wfc = w(fc). For waves, considered in this pa- 
per, when minQj.,Q* {"H} = "Hoi ^ 0. Notice, that 'H2 
in Eq. (5a) does not have aka-k and a\a*_^ terms. They 
were removed by linear canonical transformation [known 
as the Bogolubov (u, 77)-transformation] after which "^2 
takes the fully diagonal form (5a). 
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D. Three-wave interactions 



First contribution to the interaction Hamiltonian Hint 



IS 



(6a) 



fel,fe2,fe3 



describes the processes of decaying of single wave into two 
waves (1 2 processes) or confluence of two waves into a 
single one (2 1 processes). In Eq. (6) for brevity we in- 
troduced notations Oi = a^i etc. and is the Kronecker 
symbol, i.e = 1 if and only if fci + fc2 = ^a- Clearly, 



V^2- Generally speaking, also includes aia2fl3 



^2^3 

and aia2'^3 terms that describe 3 processes (conflu- 
ence of three waves or spontaneous appearance of three 
waves out of vacuum). However they can be eliminated 
by corresponding nonlinear transformation [18] that leads 
to the canonical form of "Hs, presented in Eq. (6). 

Hamiltonian H2 + %3 with Eq. (2) yields the three- 
wave equation: 



~dt 



E 

fci ,^2 



(6b) 



Two sets of terms in the RHS of this equation have 
time dependence of the form exp[— 1(^2 + ^3}t] and 
cxp[—i{LU2—LLi3)t\ correspondingly [we used shorthand no- 
tations, Wj = w{kj)]. They become important if their 
frequencies are close to the eigenfrequency of Ofe, ujk- 
LJ2 +W3 ~ or a;2 — a;3 w cj^. By relabeling the wavcvec- 
tors, we can write both of these conditions in the same 
form as follows, 



w(fci) -I- W(fe2) = ^(fes) 



(7a) 



This condition of time synchronization should be com- 
plemented by the condition of space synchronization 
that formally originates from the Kronecker symbols in 
Eq. (6b), 



fel + ^2 = ^3 



(7b) 



Both relations (7) are named the resonance conditions 
of the three- wave interactions, or conditions of the three- 
wave resonances. 

There exists a simple conditions for the three- wave res- 
onance conditions to be satisfied for the power-law dis- 
persion relations uj ^ k'^ (a=const). In 2D, it is most 
easily proved graphically, as suggested in [23]. Thus, it 
was shown that the three-wave resonance is possible if 
and only if a > 1 for the continuous case, fc e M^. Obvi- 
ously, this condition becomes a necessary condition if k is 
restricted to discrete values due to boundary conditions. 



E. Four-wave interactions 

When the three- wave resonances are forbidden, one has 
to account for processes with weaker nonlinearity, the 



four-wave interactions. The canonical part of the four- 
wave interaction Hamiltonian, 



T1I2 * * Jkl^ 
34%a2'^3a4 O34 , 



12 



(8a) 



fel,fe2,fe3,fe4 



describes a 4- wave scattering processes 2 <^ 2. Terms 
01020304 and its complex conjugate describing 4-^0 
processes can be eliminated by an appropriate nonlinear 
canonical transformation [18]. After that the four- wave 
interaction Hamiltonian takes the canonical form (8a). 
There also exist 1 <;4> 3 systems with 01020304 and its 
complex conjugate terms in V.^ [5, 6]. They can be 
treated similarly, but for simplicity we omit them in the 
present paper. 

Note that besides trivial symmetries with respect to 

the indexes permutations, 1-^2 and 3 4, the inter- 
action coefficient (8) has the symmetry Tfj = (T34) , 
because the Hamiltonian has to be real, 'Hi = %%. 

The dynamical equation for the four-wave case follows 
from (2) with the Hamiltonian "H = 7^2 + "^4: 



= WfcOfc + ^ ^ r2''3^0i0203(523 



dt 



(8b) 



fel,fe2,fe3 



Considering this equation similarly to (6b), one realizes 
that the terms in the RHS of Eq. (8b) oscillate with the 
frequencies a;2 + 1^3 — and becomes resonant if this 
combination is close to Wfc. in the other words, the con- 
dition of time synchronization (after proper renaming of 
the variables) takes the form (9a) 



oj{ki) + w(fc2) = ^(fes) -I- uj{ki) 
fei + fe2 = fe3 + fe4 , 



(9a) 
(9b) 



while Eq. (9b) represents condition of space synchroniza- 
tion that comes from the Kronecker symbol in Eq. (8b). 



F. Physical examples 

In the context of the problem at the hand, a choice 
of physically important and methodologically illustrative 
Hamiltonian systems is not an easy task. The corre- 
sponding wave systems should preferably be well-studied, 
both theoretically and experimentally (or numerically). 
They should be simple enough to be understood by the 
non-experts in the area of wave turbulence and at the 
same time not too simple in order to demonstrate the 
main characteristics of the resonant wave systems de- 
scribed by different nonlinear dispersive PDEs, with dif- 
ferent number of interacting modes and different bound- 
ary conditions. 
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Surface water waves 



3. Rossby and drift waves 



Our first example is the system of surface water waves, 
with dispersion relation of the general form: 



lgk + 



(10a) 



where g is the gravity acceleration, a is the surface ten- 
sion and p is the fluid density. For small k Eq. (10a) 
turns into dispersion law for the gravity waves: 

u^k = \fgk, (lOb) 
while for large fc it is simplified to the capillary wave form 



' a k^ 



(10c) 



In both limiting cases the dispersion law have scale- 
invariant form, ojh oc k". Notice that for the gravity 
waves a = I < 1 and therefore the leading nonlinear 
processes are four-wave scattering 2 2 with the quar- 
tets as the primary clusters, while for the capillary waves 
a = I and thus the leading nonlinear processes are three- 
wave interactions of 2 ■^4' 1 type. In this case the primary 
clusters are triads. 

Surface water waves with the general dispersion 
law (10a) can be described by the Hamiltonian equa- 
tion of motion in the canonical form (2), that turns into 
Eq. (6b) for the capillary waves and into Eq. (8b) for the 
gravity waves. Three-wave interaction coefficient for the 
capillary waves reads as 



1 _ iy/U)iUl20J3 
^23 — 



fci,— fca 



k\\/h^ ksy/k^ k2\/kik3 



where 



^/S2,/S3 = (^2 • fcs) + k2k3 . 



(11) 



The 4-wave interaction coefficient for the gravity waves 
is given by rather long expressions which can be found 
in [46]. 



2. Nonlinear Shrodinger (NLS) model 

Probably the simplest known example of the four- wave 
systems are waves in the nonlinear Shrodinger (NLS) 
model of nonlinear optical systems and Bose-Einstcin 
condensates [40, 41]. NLS waves have dispersion func- 
tion and interaction coefficient as follows: 



'■34 



1. 



(12) 



Another important example of wave system with domi- 
nating three- wave interaction, is Rossby waves, which are 
similar to drift waves in inhomogeneous plasmas. Their 
amplitudes can be described by the so-called barotropic 
vorticity equation which can be presented in the form 
similar to the canonical three- wave equation (6b), but 
all k's taking values only in half of the Fourier space, 



at 



E 



'fe2 



(fc, >0).(13) 



The phase space in this case is half of the Fourier space is 
a result is because the original equation in the a;-space is 
for a real variable (barotropic vorticity). The difference 
in the Hamiltonian structure of the Rossby and capillary 
waves yields the difference in the form of the conservation 
laws and therefore in their dynamical behavior. We will 
discuss this later in greater detail. 

Rossby waves on an infinite (or double-periodic) 13- 
plane have dispersion function [28, 47] 



I3p^k, 
1 + p2)fc2 ' 



(14a) 



where p = \/gH / f is the Rossby deformation radius H 
is the fluid layer thickness, / — 217 sin ^? is the Coriolis 
parameter, 6 is the latitude angle (/3-plane approximates 
a local region on surface of a rotating planet), is the 
planet rotation frequency and /3 is the gradient of the 
Coriolis parameter, (3 = 2Qcos9/R, and R is the radius 
of the planet. 

In the case of zero boundary conditions in a plane rect- 
angular domain [oceanic Rossby waves), the form of the 
eigenmode is given by Eq.(lc), corresponding dispersion 
function has the form 



13L 



(14b) 



Note that this dispersion relation coincides with relation 
Eq. (14a) in the limit p — >■ oo taking into account that 
{kx,ky) = (/3/2w ± Tim./L, ±Tm/L) [which follows from 
Eq.(lc)]. However, the resonant mode sets are different 
because the resonance in k is now replaced by the reso- 
nance conditions in m an n. 

One more example is Atmospheric Rossby waves, prop- 
agated on a rotating (with angular velocity fl) sphere. 

Eigenmodes in this case, I7" (sin t^. A) exp 



1) 



are proportional to the spherical functions Y^™, where 
£ < 1 and |m| < £ are integers and and A are latitude 
and longitude correspondingly. In this case dispersion 
function is of the form 



^e,m — 



2mfl 



(14c) 
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Notice that difference in the dispersion relations (14) 
leads to essential difference in the topology of resonant 
clusters, and consequently to essential difference in the 
dynamical and statistical behavior of the systems. 

For concreteness we present here the interaction coef- 
ficients of the Rossby waves in the (infinite or double- 
periodic) /3-plane [48]: 



^23 = 



Am 



k 



2y 



k 



3y 



1 + p2fc2 1 + p2fc2 1 ^ p2j^2 



The interaction coefficients for the atmospheric Rossby 
waves can be found in [39, 49, 50] and for oceanic Rossby 
waves - in [51]. 



II. REGIMES OF FINITE-SIZE WAVE 
TURBULENCE 

What happens when, due to the finite size, the num- 
ber of exact resonances and active quasi-resonances is 
depleted or absent? The finite-size effects in WT can 
be characterized by considering the nonlinear frequency 
broadening F (i.e. the inverse the characteristic time of 
nonlinear evolution) and comparing it to the frequency 
spacing A^j between the finite-box eigenmodes. For sim- 
plicity, we will restrict our attention to the periodic 
boundary conditions, in which case 



dk 



2tt 

T 



kL' 



(15) 



"Twiddle" here means that this is an order of magnitude 
relationship, which corresponds the approximate charac- 
ter of the physical estimates given below. 

The kinetic equation is applicable when F > A^^, this 
is the kinetic regime. A qualitative different behavior 
can be expected in the opposite limit F <C A^: this is 
a regime of discrete wave turbulence. These two regimes 
are realized when WT forcing is rather high (but not 
too high so that the nonlinearity is still weak) and low 
respectively. However, we will also see that there is also a 
rather wide intermediate range of forcing for which there 
is a regime with F A^j, which we will call mesoscopic 
wave turbulence. 

Name mesoscopic refers to an observation made in [35] 
and [34] that in existing numerical simulations of the 
gravity water waves there may be regimes where the sta- 
tistical properties of the infinite-box systems coexist with 
effects due to tlic^ fc-spacc discreteness associated with a 
finite computational box. In was further argued in [37] 
(in the context of MHD turbulence) that such a meso- 
scopic regime is active in a wide intermediate range of 
wave amplitudes. The key reason for such a wide meso- 
scopic range is the fact that the typical values of F for the 
discrete (dynamical) and the kinetic (statistical) regimes 
are typically strongly separated. 



A. Discrete turbulence (small box, weak waves) 

In the discrete WT regime, when F <C A^^, only the 
terms in the dynamical equations which corresponding to 
exact wavenumber and frequency resonances contribute 
to the nonlinear wave dynamics. All the other terms 
rapidly oscillate and their net long-term effect is null. 
The most clear example here is the case when there is 
no exact resonances, like in the system of the capillary 
water surface waves. In this case, the averaged (over the 
fast linear oscillations) nonlinearity is negligible and the 
turbulent cascade over scales is arrested. One can see an 
analogy with KAM theory which says that trajectories 
of a perturbed (in our case by nonlinearity) Hamilto- 
nian system remain close to the trajectories of the un- 
perturbed integrable system (in our case the linear wave 
system whose trajectories are just harmonic oscillations 
of the individual modes) if there is no resonances. Of 
course, this analogy should be taken with caution be- 
cause even in absence of the lower-order resonances (e.g. 
triad resonances for the capillary waves) higher-order res- 
onances may be important. 

Thus, for the discrete WT regime we have the following 
reduced dynamical equations 



. dak 
''~dt 



E 

1,2 



(16a) 



for the three-wave case [ Eq. (6b) in which we retain only 
exact wave resonances] and 



1,2,3 



(16b) 



for the four-wave case [ Eq. (8b) with only exact reso- 
nances left]. 

In equation (16a), factor i?^2 is equal to one when 
modes fci , k2 and fe3 are in exact wavenumber and fre- 
quency resonance, and it is zero otherwise. Respectively 
in (16b), i?34 is equal to one when modes fci, fc2, ^3 and 
ki are in exact wavenumber and frequency resonance, 
and it is zero otherwise. 

Some resonant triads/quartets (if at all present) may 
be isolated, in which case their dynamics is integrable, 
and the respective nonlinear oscillations can be expressed 
in terms of the elliptic functions. Some triads/quartets 
may be linked and form clusters of various sizes, whose 
dynamics is more complicated and to some extent may 
be chaotic, especially for larger clusters. Study and clas- 
sification of such exact resonances and their clusters was 
initiated in [39, 51, 55, 56] developed further in many pa- 
pers including [32, 34, 52-54, 57, 58]. Examples of small 
and large clusters for the Rossby waves (three-wave sys- 
tem) can be found in [52-54] and for the gravity water 
waves (four- wave system) in [58]. 

Frequency broadening F for the discrete WT can be 
estimated from the dynamical equations (16a) and (16b), 



6 



spacing A(^. Together with (15), this condition gives 



r 
r 



D 



\Wal\/^, 



(17a) 
(17b) 



\Vak\ » 



where V = V^2 s^'^d ^ = ^a^fe the interaction co- 
efScients in (16a) and (16b) respectively. Subscript D 
indicates that this is a discrete-regime estimate, and su- 
perscripts 3VF and stand for " three- wave" and " four- 
wave" respectively. Here N is the number of exact res- 
onances which are dynamically important at a fixed fc, 
which less or equal to the number of modes connected 
to k in the resonant cluster. For simplicity we assumed 
that that all the dynamically important resonances are 
local, i.e. fci ~ ^2 ^ fca ^ k. Strictly speaking, the es- 
timates (17a) and (17b) are only valid if N is not too 
large, because when M ^ 1 one should expect statisti- 
cal cancelations of the effect of different triads/quartets, 
and our estimates would have to be modified. This is the 
case for the example of MHD turbulence considered in 
[37]. Also, our estimates would have to be modified for 
systems with nonlocal in fe interactions. 

Thus, the condition of the discrete turbulence regime, 
<C A(^, becomes 



\Vak\ ^ , T Kf i for 3-wave systems, (18a) 
kLN 

\Wa\\ <^ ^\f j for 4-wave systems. (18b) 



B. Kinetic wave turbulence (infinite-box limit) 

The kinetic regime comprises the classical infinite-box 
weak WT theory, which is reviewed in the Appendix to 
this paper, including recent theory extensions to descrip- 
tion of the higher-order wave moments and PDF and 
finding solutions corresponding to turbulence intermit- 
tency [59-61]. In this regime, the frequency resonance 
broadening, denoted F^, is determined by the kinetic 
equation (AlO) for for three-wave systems and (All) for 
the four-wave systems. This gives for Tk- 

F^^"") ~ IFpnfcfcVoifc ~ \V\^\ak\^{kLf/LOk , (19a) 
F^^"") m'^nie'^/oju ^ \W\''\au\\kLr/oju (19b) 

where, for simplicity, we have assumed that the wave 
spectrum is not too narrow and the range of wavenum- 
bers interacting with k is of width ^ k. 

The upper bound for applicability of the wave kinetic 
equations follows from the condition of weak nonlinearity 
Tk -C Wfe which gives 

|yafc|(fcL)''/2<Wfc, (3-wave), (20a) 
\W\\ak\''{kLY ^ujk, (4-wave). (20b) 

Also, the wave amplitudes should be large enough for the 
broadening F^ to be much greater than the frequency 



(fcL)(''+i)/2 



YYrt > for 3-wave systems, (21a) 



^W^kV » ' fo'^ 4-wave systems. (21b) 



C. Mesoscopic turbulence and sandpile dynamics 

Consider first the case when the number of connections 
of mode k in its discrete resonant cluster is relatively 
small. A/" > 1, as it is the case, e.g., for the case of the 
gravity water waves. Comparing the range of kinetic WT, 
(21a), (21b), and the one of discrete WT, (18a), (18b), 
one can see that there exists a gap, 



km 
1 

kLN 



\Vak\ 
> > 



(fcL)(rf+i)/2 ' 

1 

{kLY+^l^' 



(3-wave), (22a) 
(4-wave). (22b) 



in which both the conditions for the kinetic WT and for 
the discrete WT are satisfied. This means that in the 
region (22a), (22b) the wave behavior is neither pure dis- 
crete nor pure kinetic WT. Existence of such a gap was 
first pointed out in [37] in the context of MHD wave 
turbulence. Region (22a), (22b) possess the features of 
both types of turbulent behavior described above. In the 
other words, in this region both types of WT may ex- 
ist and the system may oscillate; in time (or parts of the 
fe-space) between the two regimes giving rise to a qual- 
itatively new type of WT: mesoscopic wave turbulence. 
It was suggested in [34] (in the context of the surface 
gravity waves) that in forced wave systems the discrete 
and the kinetic regimes may alternate in time, see fig- 
ures 1. Namely, let us consider WT with initially very 
weak or zero intensity, so that initially WT is in the dis- 
crete regime, and let us permanently supply more wave 
energy via a weak source at small fc's. During the discrete 
phase (with fully or partially arrested cascade) the; wave 
energy accumulates until when the resonance broadening 
Tu becomes of order of the frequency spacing /\^. Af- 
ter that the turbulence cascade is released to higher fc's 
in the form of an " avalanche" characterized by predomi- 
nantly kinetic interactions. At the moment of triggering 
the avalanche, the broadening F jumps up from F = F^ 
to F = F^ ^ F/3, sec the upper figure 1. In the process 
of the avalanche release, the mean wave amplitude low- 
ers so that the value of broadening F = F^^ becomes of 
order of the frequency spacing Remember, for not 
too large N in this intermcxliate range Tk '^T^. Thus, 
at this point the system returns to the energy accumu- 
lation stage in the discrete WT regime, and the cycle 
repeats, see figures 1. Because of the obvious analogy, 
this scenario was called sandpile behavior in [34] . 

As we see, the sandpile behavior is characterized by a 
hysteresis where in the same range of amplitudes, from 
Ai to A2 in the lower figure 1, the WT intensity increases 
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FIG. 1. "Sandpilc" behavior in wave turbulence. Upper 
graph: the frequency broadening T follows the discrete tur- 
bulence dependence F = Fij until reaching the value F = A^, 
at time t = ti, at which point it jumps to the kinetic branch 
F = Tk S> Fij and rapidly drops in the kinetic regime to 
the value F = A^^ at time t = t2. Then it jumps back 
to the discrete branch F = Vd ^ F/f, after which the cy- 
cle repeats. Lower graph: the amplitude gradually grows 
to Ofe ~ Ai for t < ti and then quickly drops to A2 for 
<t <t2, after which the cycle repeats. For the three- wave 
systems Ax ~ (w/F)(ifcL)-(''+i)/2 and A2 ~ oj/{kLVM) and 
for the four-wave sy stems ~ {lj /WY''^{kL)-^'^'^+^'^'^ and 
A2 ~ ^Juj/{kLWN). 



in the discrete regime and decreases in the kinetic regime. 

For the small-amplitude part of the sandpile cycle, the 
system will be close to the critical spectrum, where res- 
onance broadening Tk is of order of the omega-spacing 
Aij. This gives the frequency spectrum cj"*', which was 
predicted in [34] and experimentally confirmed in [30] 
(c.f. for the KZ spectrum in this case [20]). Finding 
spectrum close to the large-amplitude part of the cycle 
is not so straightforward because we do not know the 
dependence oiJ\f onuj. 

So far, we only considered the case when M is not 
too large. Case A/" » 1 can be very different. Namely, 



instead of the range where both conditions satisfied si- 
multaneously, the one for the kinetic WT, (21a), (21b), 
and the one for the discrete WT, (18a), (18b), one gets 
a range where none of these two conditions are satisfied. 
This kind of mesoscopic turbulence was considered us- 
ing the MHD example in [34]. We will see that in this 
case the frequency broadening V remains of the order of 
the omega-spacing A^^ in a broad (mesoscopic) range of 
wave amplitudes. Remembering that F is a characteris- 
tic nonlinear evolution time, we note that constancy of F 
points at a possibility that the energy transfer in such a 
mesoscopic regime is driven by a hidden effectively linear 
process, which is yet to be understood. 



D. Possible coexistence of different regimes 

Strength of WT typically varies in along the turbulent 
cascade in the fc-space and, therefore, one may expect 
different wave tiirbulence regimes present in the different 
parts of the fe-space at the same instant in time. For 
example, the nonlinearity increases along the cascade to- 
ward high wavenumbers in WT of surface gravity waves 
and of MHD Alfven waves. Thus we can expect WT 
in these systems to be discrete at low fc's and kinetic 
at high fc's. Moreover, on the cross-over regions one 
can expect nontrivial gradual transition which involves 
blending and interaction of different dynamical and sta- 
tistical mechanisms. This effect is expected to be more 
pronounced if the interaction of scales is nonlocal, so that 
some wavemimber(s) from a particiilar resonant triad (or 
quartet) could be in the discrete range whereas the other 
wavenumber(s) from the same triad (or quartet) could be 
in the kinetic range. As a result, in the cross-over range 
a continuous spectrum described by the kinetic equation 
(e.g. KZ) could coexist with selected few modes belong- 
ing to isolated resonant clusters which would evolve co- 
herently at deterministic timescales. Moreover, the same 
set of modes might randomly alternate in time from be- 
ing discrete to kinetic and back, as we described above 
in the sandpile scenario. 

Some basic consequences of variability of the finite- 
size effects in the fe-space can be seen in an very simple 
kinematic cascade model suggested in [36]. This model 
builds a "cascade tree" in the following three steps. 

• Let us put some energy into a small collection of 
initial modes. We denote this initial collection of 
excited modes by S'o (e.g. in within a circle or a 
ring at small fc's which corresponds to forcing at 
large scales). One can view set as the cascade 
tree's "trunk". 

• Next, find the modes which can interact with the 
initial ones at the given level of nonlinear broaden- 
ing F. Namely, we define a new set of modes Si as 
the union of all fc's satisfying the quasi-resonance 

conditions, 

|w3 — W2 — wij < r , fe3 — fe2 — fei = (23a) 
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for the 3-wave case and 

|W4 + W3 - W2 - Wi I < r, 

fe4 + fes - fei - A;2 = 0, (23b) 

for the 4-wavc case, with all but one wavcnumbers 
in 5*0 and the remaining wavenuniber outside of 6*0. 
Provided that T is large enough, the set Si will be 
greater than ^o- Set Si comprises the cascade tree's 
"biggest branches". 

• Now iterate this procedure to generate a series of 
cascade generations Sq, Si, Sn which will mark 
the sets of active modes as the system evolves. The 
union of these sets constitutes the whole of the cas- 
cade tree with all of its bigger and smaller branches 
included. 

This model is purely kinematic. It does not say any- 
thing about how energy might be exchanged dynamically 
among the active modes, or how rapidly a certain cascade 
generation is reached. However, the kinematics alone al- 
lows one to make some interesting observations about the 
systems with variable in k finite-size effects. 

Let us consider the example of the gravity waves on 
deep water, for which the following results were obtained 
in [34]. If one starts with a set of low-fc modes, with 
broadening F below a critical value Fcrit — 1-4 x 10~^, a 
finite number of modes outside the initial region get ex- 
cited (generation 2) but there will be no quasi-resonances 
to carry energy to outer regions in further generations. 
If the broadening is larger than Fcrit, the energy cas- 
cades infinitely. Further, such the kinematic cascades 
were shown to have the fractal snowflake structure with 
the active modes being rather sparse in the front of 
the cascade propagating to higher k, with pronounced 
anisotropic and intermittent character. 

Similar picture of intermittent cascades was also ob- 
served for the capillary wave system [36]. However, be- 
cause there is no exact resonances for this system, the 
generation 1 an higher appear only if F is greater than 
some minimal value Fcrit i- Further, there exists a second 
critical value Fcrit2 > Tcriti: the number of generations is 
finite for rcrit2 > F > Fcriti and the cascade process dies 
out not reaching infinite fc's, whereas for F > Fcrit2 the 
number of generations is infinite and the cascade propa- 
gates to arbitrarily high fc's. Note that the later property 
makes the capillary wave system different from the grav- 
ity waves for which the cascade always spread through 
the wavenumber space infinitely provided F > Fcrit- 

Another example where the (three-wave) quasi- 
resonances and the kinematic energy cascades were stud- 
ied is the system of inertial waves in rotating 3D fluid 
volumes [38]. This system is anisotropic and the study 
of the kinematic cascades allows to find differences be- 
tween the 2D modes, with wavevcctors perpendicular to 
the rotation axis, and the 3D modes. It appears that the 
"catalytic" interactions which involve triads including si- 
multaneously 2D and 3D wavevectors dominate over the 
triads which involve 3D wavevectors only. 



DISCUSSION 

In this paper we have considered the three different 
regimes which can be observed in wave turbulence (WT) 
bounded by a finite box - discrete, mesoscopic and ki- 
netic. For very low amplitudes and small boxes, we 
expect the discrete WT, whose dynamics is driven by 
the exact resonances. In the opposite infinite-box limit, 
we expect the kinetic WT, which is driven by quasi- 
resonances and for which the exact resonances do not 
play a role as they are hugely outnumbered by the quasi- 
resonances. This is the classical and the most studied 
WT regime, and it is summarized in our Appendix. In 
the middle, there is a regime of the mesoscopic WT. We 
have shown that this regime is characterized by sandpilc- 
like oscillations between the discrete and the kinematic 
regimes (if the size of the active resonant clusters is small) 
or it settles to an intermediate (critical) state in which 
the nonlinear frequency broadening is of order of the fre- 
quency spacing between the discrete modes for a wide 
range of wavenumbers (if the size of the active resonant 
clusters is large). 

The key fact that has led us to the observation that 
the mesoscopic regime should be realized in a wide range 
of wave intensities, is that the dependence of the fre- 
quency broadening on the wave intensity is very different 
for the dynamical and the kinetic equations; cf Fj^ given 
by (17a), (17b) and Tk given by (19a), (19b). Thus, 
for the same wave intensities in which F/j and Tk are 
typically very different in size, and there exist a wide 
mesoscopic range where either both the discrete and the 
kinetic regimes can exist, or none of them is realizable, - 
hence the two types of the mesoscopic behavior described 
above. 

Signs of bursty behavior typical of the sandpile behav- 
ior suggested in this paper has already seen in laboratory 
and numerical experiments [30, 34]. In future, one should 
aim to perform more direct diagnostics of the quantities 
allowing to identify and to distinguish the different WT 
regimes describe in the present paper, including the non- 
linear frequency broadening and character of its evolution 
in time. 



Appendix A: Kinetic wave turbulence 

Classical Wave Turbulence theory provides a statisti- 
cal description of weakly nonlinear waves with random 
phases. As discussed above, theory of wave turbulence is 
valid in a range of wave-field strengths such that 

1 > - > — ~ T^, (Al) 

u)k tOk kL 

where F is given by (19a) or (19b) for the three- or four- 
wave processes respectively. 

The most popular statistical object in the theory of 
wave turbulence is the waveaction spectrum, although 
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theory of wave turbulence has been recently extended to 
description of higher moments and probability density 
functions (PDF) in [59-61]. This allowed to deal with 
non-gaussian wave fields, as well as to study validity of 
the underlying statistical assumptions such as e.g. ran- 
dom phases. We will now briefly describe these results. 

Let us represent the complex amplitudes as = 
\/jfc7/'fe with wave intensity Jk € M"*" (positive real num- 
ber) and phase factor V'/c € (complex number of length 
1). Let us define the M-mode joint PDF p(^) so that 
the probability for the wave intensities of the selected M 
modes, Jk, to be in the range (s^, Sk + dsk) and for their 
phase factors if^k to be on the unit-circle segment between 
a and ^k+d^k is P^^^ IlfeLi dsk \d^k\- (Therefore p(^) 
is a function of 2M + 1 variables: M amplitudes, M 
phases and time). 

Notion of random phases refers to the cases where all 
factor V'fc are statistically independent and uniformly dis- 
tributed on i.e. 



-p{M) 



(2^) 



-p(M) 



M a 



(A2) 



for any M < N, where N is the total number of dynam- 
ically active modes. Here vi^^ is the joint PDF of the 
amplitudes only. Kinetic WT considers wavefields with 
random phases at some initial time and with intensities 
satisfying condition (Al). This leads to the following 
equation for the joint PDF for the three-wave case. 



g-piN) 

dt 

X 



5_ 

5s 



167r j |V23p(J(wi -i02- W3)(5(fci - k2- ks) 
p(^) ) dkidk2dk3, (A3) 



S1S2S3 



Ss 



where 



-. This equation was first 



5] _ _6 S__ 

6s 5si 5s2 Sss 
derived for a specific example of waves in anharmonic 
crystals by Peierls [19] and for general three- wave sys- 
tems in [59, 61, 62]. It was also extended to the four wave 
systems in [60] . Note that the phase variables are not in- 
volved in these equations. Therefore, the random phase 
assumption is consistent with these equations, namely 
the system which has random phases initially will remain 
random-phased over the typical nonlinear time (i.e. its 
PDF will remain independent of ^'s). Thus, equations 
for the joint PDF (A3) allows an a posteriori justifica- 
tion of the random phase assumption underlying their 
derivations. 

I 



However, as we already mentioned, the most frequently 
considered object in the theory of wave turbulence is the 
spectrum which is defined as 



nk 



(A4) 



where d is the dimension of the space and the angu- 
lar brackets mean the ensemble averaging over the wave 
statistics. The spectrum is a one-mode statistical object, 
and it is the first in the series of one-mode moments. 



27r 

T 



)pd „oc 



Sk)dsk. 



Note that for deriving closures for the one-mode ob- 
jects the random phase property is insufficient and one 
has to assume additionally that the amplitudes Jk are 
also statistically independent of each other at diff'erent 
fc's. Statistical independent of the amplitude can also be 
justified based on the equation for the joint PDF (A3), 
although this issue is more subtle than the phase ran- 
domness because variables Sk do not separate in (A3) 
and, therefore any product factorization of the joint PDF 
in terms of the one-mode PDF's would not generally be 
preserved by the nonlinear evolution. However, this sit- 
uation seems to be typical for many systems, e.g. for 
the relation between the multi-particle and one-particle 
distribution functions described by the Louisville and 
Boltzmann equations respectively. In these situations, 
a sufficient for the closures property is that the low- 
order PDF's, P(^) with M <^ N, can be product fac- 
torized. It can be seen from (A3) that it is the case for 
the weakly nonlinear wave systems, i.e. that factoriza- 
tion, V^'^'> = U^^i'Pk ^ + 0{M/N), survives over the 
characteristic nonlinear time. 

Importantly, the distribution of wavefields in the the- 
ory of kinetic WT docs not need to be Gaussian or 
close to Gaussian, and one can consider evolution of 
the one-mode PDF's V^^^ that correspond to strongly 
non-gaussian fields (Gaussian fields would mean V^^^ <^ 
e~^l ^■^^). Integrating the joint PDF equation (A3) we get 



dt dsk ' 
with F is a probability fiux in the s-space, 



Fk = -sk{,kri'^ + nk^). 



5sk 

where for the three-wave case we have: 



(A5) 



(A6) 



% = 47r j {\Vi2\'^S{cok - wi - 0J2)S{k - fei - fe2) + 2\Vki\'^S{L02 - ojk - 0Ji)6{k2 -k- fei))nin2 dkidk2, (A7) 
7fe = Stt j {\Vi2f5{iOk - wi - cJ2)S{k - fei - fe2)ni -|- \Vki\^Ali5{cj2 - Wfe - wi)(5(fe2 - fe - fei)(ni - 712)) dkidk2. 
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Equation (A5) has an obvious exponential solution which 
corresponds to a zero flux F: 



>(i) 



1 



(Jk) 



-Sk/(Jk) 



which corresponds to Gaussian statistics of the wave field 
Ofc. However, there are also solutions corresponding to 
F =const7^ which for Sk ^ {Jk) has a power-law 
asymptotic [59, 60], 

>(i) - 



K' = - 



IkSk 



These solution corresponds to enhanced probability (with 
respect to gaussian) of strong waves which is called in- 
temittency of WT. Here, the constant flux in the am- 
plitude space F can be associated with a wavebreaking 
process the exact form of which depends on the physi- 
cal system. For example, for the gravity water surface 
waves the wavebreaking process takes form of whitecap- 
ping, and for the focusing NLS system the wavebreaking 
is represented by filamentation/coUapsing events. Ob- 
viously, this power-law tail of the PDF cannot extend 



to infinity because the integral of the PDF must con- 
verge. Thus, there exists a cutoff which can also be as- 
sociated with the wave breaking, which can simply be 
understood that the probability of waves with amplitude 
greater than some critical value must be zero. Such criti- 
cal value roughly corresponds to the amplitude for which 
the nonlinear term becomes of the order of the nonlinear 
one so that the WT description breaks. 

Multiplying equation (A5) by sJJ and integrating over 
Sk, we have the following equation for the moments 



M. 



(p) 



dt 



(p-i) 



(A8) 



which, for p = 1 gives the kinetic equation for the wave- 
action spectrum, 



di 



nk = -Iknk + Vk- 



(A9) 



Substituting into this equation expressions for and rjk , 
we obtain more familiar forms of the kinetic equations: 



di 



rifc = 47r j |142p<5('^fc - oji- co2)S{k - fei - A;2)(nin2 - niUk - n^rik) dk^dk^ 
-HStt j \ V^i\'^5{ijj2 - Wfe - LOi)5{k2 - k - ki){nin2 - fiiUk + n2nk) dkidk2, 



and for the four- wave case 
d 



dt 



nk 



23 +u!i-u!2- uj3)nknin2n3S{k + ki-k2- fea) 



Vrifc 111 712 "3 



(AlO) 



dkidk2dki . (All) 



Based on Eqs. (AlO) and (All) [or Eqs. (A5) and (A8)] 
one can obtain the estimate for the nonlinear frequency 
broadening in the WT regime, i.e. inverse characteristic 
time of the nonlinear evolution as in (19a) or (19b). This 
leads to the WT applicability condition (21a) or (21b). 

Classical statistical approach allows to obtain some 
interesting and physically relevant solutions, such as 
Kolmogorov-Zakharov (KZ) spectra corresponding to the 
energy and waveaction cascades through scales. Such 
solutions can be obtained analytically using so-called 
Kraichnan-Zakharov transformation, as well as from the 
scalings of the frequency and the interaction coefficients 
based on the dimensional analysis. Discussion of these 
issues is beyond the scope of our review, and the inter- 
ested reader is referred for details to book [18]. Here, it 
suffices to say that in most systems there exists a short- 
cut way to obtain KZ spectra. It works for the systems 
with only one relevant dimensional parameter, for exam- 
ple the gravity constant g for the water surface gravity 
waves, surface tension constant a for the capillary waves. 



speed of sound Cg for acoustic turbulence, quantum of 
circulation k for Kelvin waves on quantized vortex lines, 
etc. In this case the ID energy spectrum Ek ~ k'^ can 
be immediately obtained from the physical dimension of 
this constant which gives for the direct cascade [36]: 

, ^ 5-3Q!-d , 
t/ = 2aH-d-6-h , (A12) 

where a is the power of the dispersion relation uj ^ k" 
(which is uniquely determined by the above dimensional 
constant), d is the dimension of the system and N is the 
number of waves involved in the resonance interaction. 
For example, for the water surface gravity waves we have 
Ek ^ for the capillary waves Ek ~ k~'^^^ (both 

of these spectra are called Zakharov-Filonenko spectra 
[20, 21]), for acoustic turbulence Ek ~ (Zakharov- 
Sagdecv spectrum [29]). For Kelvin waves on quantized 
vortex lines, considering them as a local six-wave process, 
one formally gets Ek ^ fc"^/^ (Kozik-Svistunov spectrum 
[3]). However, this spectrum was recently shown in to be 
nonlocal 
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Similar approach one can use for finding the inverse cascade spectra, e.g. for the water surface gravity waves 

or Kelvin waves [36]. 
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